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1.  INTRODUCTION  ^ 

The  accurate  evaluation  of  the  dynamic  fracture  toughness  of  advanced  materials  is 
becoming  increasingly  important.  Under  dynamic  fracture,  material  at  the  crack  tip  is 
strained  suddenly,  and  if  is  rate  sensitive,  may  offer  more  resistance  to  fracture  than  at  quasi¬ 
static  strain  rates.  High  stain  rate  fracture  testing  is  of  interest  because  many  structural 
components  are  subject  to  high  loading  rates  in  service  or  must  survive  high  loading  rates 
during  accident  conditions.  Thus,  these  components  must  be  designed  against  crack 
initiation  under  high  loading  rates  or  to  arrest  a  rapidly  propagating  crack.  The  fracture 
toughness  of  a  material  loaded  suddenly  is  generally  higher  than  when  the  load  is  applied 
quasi-statically;  therefore,  dynamic  toughness  is  an  effective  design  parameters. 

While  there  exists  standardized  method  for  determining  the  fracture  toughness 
parameters  Kj^  and  under  quasi-static  loading,  to  date,  a  umfied  view  on  the  best 
measurement  procedure  to  determine  ihe  dynamic  fracture  toughness  of  advanced  materials 
has  not  been  achieved.  Hence,  researchers  still  need  to  identify  a  simple  and  reliable 
experimental  method  capable  of  providing  accurate  dynamic  fracture  toughness 
measurements.  Further  understanding  of  the  failure  phenomena  (such  as  crack  propagation, 
material  instabilities,  etc.)  in  advanced  layered  and  nano-materials  require  the  development  of 
novel  experimental  techniques.  The  motivation  of  the  present  investigation  is  to  establish  a 
reliable  experimental/computational  procedure  to  obtain  both  accurate  dynamic  fracture 
toughness  and  crack  dynamic  propagation  parameters. 

We  are  investigating  failure  modes  by  means  of  the  impact  of  pre-cracked  4-point 
bending  specimens.  Two  fidl  field  surface  displacement  measurements,  dynamic  white-light- 
speckle  and  coherent  gradient  sensing,  are  being  used  to  examine  strain  fields  at  crack-tips  in 
ceramics  and  layered-materials.  A  light  gas  gun  and  a  Kolsky  bar  are  employed  to 
dynamically  load  the  samples.  Real  time  photographs  are  recorded  using  the  Gordin  Model 
220  high-speed,  8-CCD,  camera  purchased  under  this  grant.  In  combination  with  the 
measurement  of  dynamic  load  history,  crack  propagation  using  a  crack  gage,  the  following 
fracture  parameters  can  be  obtained  from  these  full  field  tip  images:  1)  dynamic  fracture 
toughness  at  initiation,  2)  crack  propagation  speed  and  3)  stress  intensity  factor  history. 


2.  RESEARCH  METHODOLOGY 


A-stored-energy  Kolsky  bar  as  shown  Figure  1  is  employed  to  establish  and  calibrate 
the  dynamic  firacture  measurement  technique.  Figure  2  shows  the  Single-Edge  Notched 
Beam  (SENB)  samples  used.  During  the  test  process,  the  SENS  samples  are  mounted  in  the 
Kolsky  bar  and  pre-loaded  to  keep  them  on  position.  The  5  Watt  Verdi  laser  (Coherent  Co.) 
is  used  to  illuminate  the  sample  using  an  optical  beam  splitter  and  a  light  diffuser.  A  Model 
K2  Infinity  Long-Distance  Microscope  and  the  newly  purchased  high  speed  8-CCD  Cordin 
Camera  model  220  (  as  shown  in  Figure  3),  with  up  to  100  million  frames  per  second,  images 
the  crack  tip  during  sample  deformation  and  failure.  The  high  speed  camera  is  triggered  by  a 
Tektronix  model  TDS520C  oscilloscope  with  a  frequency  response  up  to  500MHz  and  a 
maximum  sampling  rate  of  IGs/s.  The  incident  and  transmitted  loads  are  measured  by  two 
high  frequency  response  load  ceUs.  The  load  cells  are  Quartz  Load  Washer  Type  9001  from 
Kisder  Co.  with  a  range  of  35KN.  The  load  cells  are  used  wifii  a  pre-load  of  10%  the 
maximum  load  to  ensure  stability  and  rapid  response.  A  Charge  Amplifier  Model  501  OB 
from  Kistler  Co.  with  a  cut  off  frequency  of  180  KHz  is  also  utilized  to  energize  the  load 
cells.  The  crack  initiation  time  is  measure  by  a  crack  gauge  that  is  glued  on  the  sample.  A 
schematic  of  the  experimental  setup  is  shown  in  Figure  4. 

For  the  full  field  measurements,  two  fid!  field  measurement  techniques  are  now 
being  used,  Speckle-Direct  Image  Correlation  and  coherent  gradient  sensing  (CGS) 
techniques.  When  the  crack  tip  display  little  plastic  deformation,  the  Speckle-DIC  method  is 
ideal.  By  contrast,  when  large  plasticity  develops  at  the  crack  tip,  CGS  is  used. 


Figure  1:  Picture  of  the  stored-torque  torsional  Kolsky  Bar  Apparatus 


Figure  2:  Single-Edge  Notched  Beam  (SENB) 


Figure  3:  Gordin  Model  220  high-speed  8-CCD  camera 
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Figure  4:  Schematic  description  of  the  experiment  setup 


2.1  Speckle-DIC  Method 

Figure  5  shows  the  detail  Speckle-DIC  setup  in  our  Dynamic  Inelasticity  Laboratory. 
Before  the  test,  the  sample  is  mirror-polished.  A  random  black-and-white  speckle  pattern  is 
then  produced  on  one  surface  of  the  sample  and  the  crack  gage  is  glued  on  the  opposite 

Several  tests  were  performed.  Figure  6  shows  the  typical  speckle  images  of  SENB 
Aluminum  specimens  obtained  using  die  present  setup.  The  digitized  images  are  currendy 
being  analyzed  with  digital  correlation  procedures  to  obtain  the  displacement  and  strain 
fields.  With  this  full  field  information,  in  combination  with  crack  initiation  time,  measured 
by  the  crack  gage,  and  FEM  simulations,  the  d3mamic  firacture  toughness  and  stress  intensity 
factor  history  of  interest  during  the  whole  firacture  process,  will  be  obtained.  We  expect  to 
correlate  several  experimental  measurements  in  order  to  constrain  the  numerical  simulations 
that  are  used  to  estimate  K„)  or  J  .  In  the  case  of  britde  materials,  the  crack  tip  position 
history  can  be  clearly  seen  in  the  images.  From  the  tip  positions,  the  average  crack 
propagation  speed  can  be  estimated  by  dividing  the  relative  distance  of  two  near  crack  tip  by 
the  time  difference. 


Figure  5;  Speckle  set-up 


Picture  3:  300ms 


Picture  4:  500ms 


Picture  5:  700ms 


Picture  6:  900ms 


Figure  6:  Speckle  images  of  SEND  Aluminum  6061-T6  Specimen 


2.2  Coherent  Gradient  Sensing  Method  (  CGS) 

For  some  material large  amounts  of  plasticity  develop  near  the  crack  tip.  CGS  is  the 
suitable  full  field  measurement  method  in  diis  case.  Figure  7  shows  the  schematic  setup  of 
our  present  CGS  test.  The  5  Watt  Verdi  laser  (Coherent  Co.),  optical  components,  specimen, 
and  the  Cordin  Model  220  high  speed  camera  are  aligned  precisely.  The  laser  light  must  be 
perfecdy  perpendicular  to  the  sample  to  obtain  a  working  axis  over  which  the  optical  setup  is 
mounted.  All  CGS  optical  components  are  layed  on  that  optical  path.  The  sample  surface  is 
mirror-polished  to  achieve  the  desired  reflectivity.  The  crack  gage  is  glued  on  the  back 
surface  of  the  sample  to  provide  a  redundant  measurement  of  crack  speed. 

Several  tests  were  performed.  Figure  8  shows  typical  high  speed  images  of  SENB 
PMMA  specimens  during  the  dynamic  crack  propagation  process.  From  the  images  recorded 
with  the  high  speed  camera,  we  are  now  obtaining  in-plane  gradients  of  out-of-plane  surface 
displacements  around  the  crack  tip.  Since  it  is  relatively  insensitive  to  vibrations  and  rigid 
body  motions,  this  technique  has  potential  for  dynamic  crack  growth  measurements  in 
advanced  materials. 


Figure  7:  Schematic  CGS  Setup 


Picture  3:  200ms  Picture  4:  300ms 


Figure  8:  High  Speed  Camera  Images  of  PMMA  SEND  Specimen 


3.  DATA  ANALYSIS 


When  a  SENB  sample  is  rapidly  loaded,  diere  are  three  typical  time  domains:  short, 
intermediate,  and  long.  In  the  short  time  domain,  there  is  a  complex  wave  propagation  in  die 
specimen  and  inertia  effects  dominates  the  specimen  response.  In  this  short  time,  reflecting 
stress  waves  that  pass  dirough  the  specimen  constructively  and  destructively  interfere  with 
one  anotiier,  resulting  in  a  highly  complex  time-dependent  stress  distribution.  It  is  not  clear 
whether  or  not  K  or  J-dominated  field  exists  at  this  stage.  In  the  intermediate  time  domain, 
the  global  inertia  effects  are  significant  but  local  oscillations  at  the  crack  are  small.  In  the 
long  time  domain,  the  deformation  dominates  the  specimen  and  K  or  J-dominated  fields 
should  exist.  Nakamura  et  al.  (1986)  defined  a  transiHon  time,  ,  when  the  kinetic  energy  and 
the  deformation  energy  are  equal.  Inertia  effects  dominate  prior  to  the  transition  time,  but 
the  deformation  energy  dominates  at  time  significantly  greater  than  .  In  the  latter  case,  a  K 
or  J-dominated  field  should  exist  near  the  crack  tip  and  quasi-static  relationships  can  be  used 
to  infer  K  or  J  from  global  load  and  displacement  measurements.  However,  during  the 
transition  time,  quasi-static  relationships  can  not  be  used  if  crack  propagation  occurs  before 
.  The  wave  propagation  phenomenon  is  a  function  of  the  material  ductility.  Brittle  and 
fiber  composite  materials  may  not  achieve  a  K  or  J-dominated  condition.  Hence,  in-depth 
investigations  are  required  to  gain  insight  into  these  features. 

Based  on  elastic  BemouUi-Euler  beam  theory,  Nakamura  et  al.(1986)  proposed  a 
formula  for  the  transition  time,  viz.. 


Where  D  is  a  displacement  coefficient  defined  by 


D  = 


A(0 


A  is  a  geometry  factor ,  W  is  the  specimen  width,  Cois  the  longitudinal  wave  speed  of  the 


materials,  A  is  the  loading  line  displacement,  and  ^  (0  is  the  displacement  rate.  For  the  3- 
point  bending  specimen,  A  is  given  by 


A  = 


SEBC 

W 


If  fracture  occurs  after  2  ,  then  it  is  normally  thought  that  the  static  formulas  for 

fracture  toughness  can  provide  an  accurate  quantification  of  J,c  or  Kjc. 

It  is  worth  mentioning  that  for  ductile  materials,  the  transition  time  requirement  can 
be  met  by  decreasing  the  displacement  rate  or  the  width  of  the  specimen,  but  for  brittle 
materials,  this  requirement  is  difficult  to  meet.  It  is  obvious  that  the  static  formula  can  not  be 
used  if  the  d3mamic  fracture  properties  of  running  cracks  and  the  effect  of  high  loading  rated 
on  fracture  are  of  interest.  In  those  circumstances,  a  dynamic  data  interpretation  method 
should  be  used  to  de-convolute  experimental  test  data.  In  the  experimental  process,  there  are 


significant  differences  between  ductile  and  britde  materials.  So  we  will  discuss  them 
separately. 

3.1  Ductile  materials 

For  ductile  materials,  it  has  been  su^ested  that,  provided  a  one  parameter 
representation  of  the  crack  tip  fields  remains  valid,  dynamic  firacture  toughness  can  be  a 
condition  for  onset  of  crack  growtli.  If  the  plasticity  near  the  crack  tip  is  small,  the  Kp  can 
be  an  alternative  of  J^.  To  date  relatively  little  experimental  work  has  been  done  on 
determining  fracture  parameters  for  ductile  fracture  under  dynamic  loading  conditions. 
Although  some  attempts  have  been  made,  the  data  inteipretation  is  still  based  on  quasi-static 
formulae  for  J.  Simple  and  reliable  test  methods  for  characterizing  dynamic  firacture 
toughness  in  metals  has  not  been  established. 

Theoretically,  the  transition  time  requirement  can  be  met  by  decreasing  the 
displacement  rate  or  the  width  of  die  specimen.  It  is  obvious  that  this  requirement  can  not 
be  met  if  the  dynamic  fracture  properties  of  fast  propagating  cracks  or  high  loading  rates  are 
of  interest.  In  such  cases,  dynamic  data  interpretation  methods  should  be  used.  In  this 
situation,  the  dynamic  firacture  toughness  can  be  estimated  by  relating  it  to  the  measured 
crack  tip  opening  displacement  (CTOD)  history.  The  details  of  die  analysis  procedure  are 
the  following: 

a.)  Asymptotic  elastic-plastic  crack  tip  field 


Hutchinson  (1968)  and  Rice  and  Rosengren  (1968),  collectively  referred  to  as  HRR, 
considered  the  case  of  monotonically  loaded  stationary  crack  in  a  material  described  by  a  J^- 
deformation  theory  of  plasticity  and  a  power  law  hardening  relationship  between  the  plastic 
strain  and  stress.  They  showed  that  the  strain  components  in  the  crack  tip  region  scale  with 
the  value  of  the  J-integral.  Within  a  small  strain  assumption,  the  asymptotic  solution  of  the 
elastic-plastic  field  equations,  in  the  crack  tip  region,  has  the  form 


tj  o 


-in/(n+l) 


Eijin,o) 


<^ij 


nn/(rt+l) 
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as  r— »0.  CT(,is  the  tensile  yield  stress,  fi^is  the  equivalent  tensile  yield  strain,  n  is  the 
hardening  exponent.  The  dimensionless  quantity, /„,  is  defined  by  Hutchinson  (1968).  The 

amplitude  factor  J  is  the  value  of  Rice's  /-integral  (Rice,  1968).  It  has  been  su^psted  that, 
provided  a  one  parameter  representation  of  the  crack  tip  fields  remains  valid,  a  condition  for 
onset  of  crack  growth  is  the  attainment  of  a  critical  value  of  /. 

b.)  Crack  tip  opening  displacement  (CTOD)  associated  with  the  HRR  singular  field 


If  the  CTOD  is  defined  using  the  intersection  of  a  90"  vertex  with  the  crack  flanks 
as  shown  in  Figure  9  ( Shih,  1981),  the  J-integral  can  be  related  to  CTOD  using  the  HRR 
singularity  field  as 


Crack  face 


7  = 


— 5 
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Figure  9;  CTOD  defined  by  Shih  (1981) 


where  ,  5  is  the  CTOD,  J  is  the  value  of  J-integral,  CT^is  the  yield  stress  and  d„  is  a  material 
dependent  dimensionless  constant  as  defined  by  Shih.  (1981).  For  dynamic  loading,  the 
above  equation  becomes 


7(0  = 


<^o 


Sit). 


At  crack  initiation  time,  which  is  defined  as  the  crack  propagation  time. 

The  CTOD  can  be  obtained  either  experimentaUy  or  numerically. 

The  first  method  is  to  obtain  the  CTOD  at  crack  initiation  time  directly  from  high 
speed  camera  images.  Typical  speckle  pictures  are  shown  in  Figure  5.  If  at  the  crack  initiation 
time  the  CTOD  is  measured,  the  can  be  obtained  direedy  firom  the  above  formula. 

Another  method  to  identify  CTOD  is  using  FEM  simulations  with  the  incident 
loading  and  transmitted  loading,  recorded  with  load  cells,  as  input  loading  conditions.  The 
CTOD  defined  above  can  be  obtained  from  the  FEM  solutions  and  therefore  the  can 
then  be  obtained  with  the  numerically  CTOD  at  the  crack  initiation  time  measured 
experimentally.  Alternatively,  a  calculation  of  J  can  be  performed  and  7^  =  7  {t- )  computed. 


3.2  Brittle  materials: 

For  brittle  materials,  it  is  very  difficult  to  meet  the  transition  time  requirement. 
Hence,  the  dynamic  data  analysis  should  be  performed  in  almost  all  the  time  domain. 
Plasticity  is  supposed  to  be  very  small  close  to  the  crack  tip  and  a  dynamic  stress  intensity 
factor  K  is  supposed  to  control  the  crack  behavior.  In  the  present  study,  the  dynamic 


fracture  toughness  will  be  estimated  by  correlating  the  displacement  from  experiment 
with  those  from  asymptotic  solution. 

a.)  Asymptotic  elastic-plastic  crack  tip  field 

The  asymptotic  crack  tip  field  under  steady-state  conditions  is  given  as  an  infinite 
series  by  Nishiol^  et  al.(1983).  The  displacement  fields  under  Mode  I  condition  can  be 
written  as 


U  =  - -J ^  -  h{n)ri  cos(-^)  } 
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C  is  the  speed  of  the  propagating  crack  tip,  K"  =  (3— V)/(l-hV)  for  the  plane  stress 
conditions,  V  is  Possion's  ratio,  p  is  the  shear  modulus  of  elasticity  and  p  is  the  density.  In 
the  series  solution,  the  coefficient  of  the  first  term  K,  is  equal  to  K(t).  At  the  crack  initiation 
time,  ATo  = 

A  region  on  the  specimen  in  which  the  most  singular  term  of  the  asymptotic 
expansion,  for  the  above  displacements,  describes  the  actual  displacement  field  is  called  a  K- 
dominant  region.  If  experimental  displacement  data  are  taken  at  point  (x,  y)  exhibiting  such  a 
K-dominant  region  then  extraction  of  Kp  can  be  obtained  by  performing  a  least-squares 
fitting  procedure  of  the  above  equation.  FEM  simulation  can  be  employed  to  decide  the  near 
crack  tip  three-dimensional  region.  When  analyzing  the  experimental  displacement  data,  data 
points  lying  within  this  region,  are  excluded  since  this  is  a  three-dimensional  deformation 
zone. 


Although  the  above  asymptotic  field  has  the  steady-state  assumptions  (i.e.,  constant 
crack  tip  speed),  the  assumption  is  acceptable  in  deciding  the  dynamic  fracture  initiation 

toughness. 

b.)  Determination  of  displacement  fields  experimentally 

In  order  to  use  the  above  solutions  to  identify  the  dynamic  fracture  toughness, 
displacement  fields  must  be  measured.  At  present,  two  methods  will  be  used:  speckle-DIC 

and  FEM  simulation.  ,  •  -n  l 

In  the  speckle  method,  a  subset,  such  as  10X10  mm ,  near  the  crack  tip  will  be 

chosen  and  will  be  correlated  digitally  from  one  image  to  another.  Once  the  location  of  this 
subset  in  the  deformed  images  is  found,  displacement  of  this  subset  is  known.  The  dynamic 
fracture  toughness  is  obtained  through  comparing  the  DIG  displacement  vtith  asymptotic 
displacement  field.  The  fracture  toughness  is  then  obtained  at  the  crack  initiation  time. 

In  the  FEM  simulation  method,  the  incident  load  and  the  transmitted  load  obtained 
from  the  experiment  are  used  as  the  FEM  input  load  conditions.  FEM  simulations  are  then 
used  to  obtain  the  displacement  field.  The  Dynamic  fracture  toughness  is  obtained  through 
comparison  of  the  FEM  displacement  with  measured  displacement  field  at  initiation  time. 
The  fracture  toughness  is  then  obtained  at  the  crack  initiation  time. 

The  present  metiiodology  also  allows  the  study  of  running  cracks  at  constant  speed 
(steady  state)  or  even  tiirough  FEM  simulations  the  Kq  quantification  for  the  case  of 
accelerating  cracks. 

4.  FUTURE  WORK 

The  present  method  will  be  further  developed  so  that  it  can  be  employed  to  study  dynamic 
failure  mechanisms  of  nano-materials  and  layered  materials  under  various  dynamic  loading 
conditions  and  various  environmental  conditions  (high  and  low  temperatures).  The  effect  of 
material  microstructure  on  the  dynamic  failure  of  advanced  materials  will  be  studied  using 
the  methodology  previously  discussed. 

New  techniques  Aat  can  trace  the  development  and  propagation  of  material 
instabilities  in  advanced  materials  will  also  be  developed  combining  high  speed  photography, 
laser  and  optical  techniques.  This  will  help  the  understanding  of  mechanisms  leading  to  shear 
band  initiation,  propagation  and  related  failure  phenomena.  With  these  investigations,  shear 
band  initiation  and  propagation  criterion  of  advanced  materials  will  be  established. 

5.  PERSONNEL 

Principal  Investigaton  H.D.  Espinosa 

Graduate  Students:  H.  V.  Arrieta,  A.  Patanella,  H.  Zhang 
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